Absence of Reentrance in the Two-Dimensional XY-Model with Random Phase Shift 



in 



Thomas Nattermann^''^ Stefan Scheidl^''^ Sergey E. Korshunov'^^ and Mai Suan 
Institut fiir Theoretische Physik, Universitdt zu Koln, D-50937 Koln, Germany 
Landau Institute for Theoretical Physics, Kosygina 2, 117940 Moscow, Russia 
Institute of Physics, 02-668 Warsaw, Poland 
(February 1995) 

We show, that the 2D XY-model with random phase shifts exhibits for low temperature and 
small disorder a phase with quasi-long-range order, and that the transition to the disordered phase 
is not reentrant. These results are obtained by heuristic arguments, an analytical renormalization 
group calculation, and a numerical Migdal-Kadanoff renormalization group treatment. Previous 
predictions of reentrance are found to fail due to an overestimation of the vortex pair density as a 
consequence of independent dipole approximations. At positions, where vortex pairs are energetically 
favored by disorder, their statistics becomes effectively fermionic. The results may have implications 
for a large number of related models. 
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We reconsider in this paper the 2-dimensional XY- 
model 



n = -j 



<ij> 



COS 



Aij) 



(1) 



with quenched random phase shifts Aij on the bonds, 
where i,j run over the sites of a square lattice. For sim- 
pHcity we assume, that the Aij on different bonds are un- 
correlated and gaussian distributed with mean zero and 
variance a. 

Model (1) describes e.g. 2-dimensional XY-magnets 
with random Dzyaloshinskii-Moriya interaction [Q. 
Other realizations are given by Josephson-junction ar- 
rays with positional disorder |^ and model vortex glasses 
In particular, in the case of the so-called gauge glass 
model, one assumes Aij to be uniformly distributed be- 
tween and 27r. We expect, that our model with gaus- 
sian disorder is equivalent to the gauge glass model when 
a oo. 

For vanishing Aij model (1) undergoes a Kosterlitz- 
Thouless (KT) transition, at which the spin-spin corre- 
lation exponent 77 jumps from 1/4 to zero 

Weak disorder, tr ^ 1, should not change much this 
picture. In the spin wave approximation one obtains rj 
— -^{T / J + cr), which remains now finite even at T = 0. 
The features of the KT-transition are essentially pre- 
served, but the transition is shifted to lower tempera- 
tures and the jump of 77 at the transition is diminished 
iQ. The actual transition temperature Tc{a) < 7+(cr) 
depends on the bare value for the vortex core energy Ec, 
here T± = f J[l±(l - Sct/tt)^/^]. In the limit 00, 
Tc = T+. 

Strong disorder will suppress the quasi-long-range or- 
der of the KT phase [|. In particular, ii Q — 
1^ S<p/a(j> ^ij order one, vortices are gener- 



<plaq> 



denotes 



ated even at zero temperature. Here ^ 
the sum over the four bonds of an elementary plaquette. 



Rubinstein, Shraiman and Nelson (RSN) Q extended 
the Coulomb gas description of the KT-transition Q to 
the presence of randomly frozen dipoles arising from the 
random phase shifts. Surprisingly, they found a second 
(reentrant) transition at Tre{<y) (< r_(tT)) to a disor- 
dered phase at low temperatures (see Fig. 1). Tre{<j) 
bends towards higher temperatures for increasing disor- 
der. The two lines T± merge at CTc = 7r/8. For a > <Tc 
there is no ordered phase. The precise value of rre(cr) 
depends again on Ec- Similar results were obtained in 
Ref. [2]. 

Korshunov [5] has argued, that the intermediate phase 
in the range Treic) < T < Tc{a) with quasi- long-range 
order is probably not stable, if in addition to the screening 
of Coulomb charges by neutral pairs of charges, consid- 
ered in [1], screening by larger complexes of charges in 
different replicas are taken into account. 

Experiments [6] as well as Monte Carlo studies [7] in- 
dicate no reentrance. Also, Ozeki and Nishimori [8] have 
shown for a general class of random spin systems, which 
include (1), that the phase boundary between the KT- 
and the paramagnetic phase is parallel to the tempera- 
ture axis for low T. Thus they exclude a reentrant tran- 
sition, provided the intermediate KT phase exists. How- 
ever, they cannot rule out the possibility, that the KT- 
phase disappears completely, as suggested in [5]. 

We will argue below, that the reentrant transition is 
indeed an artefact of the calculation scheme used in [1], 
[2] and that the KT-phase is stable at low temperatures 
with Tc{a) ^ for (T ^ tt/S (see Fig.l). Since the 
renormalization group (RG) flow equations (4) (see be- 
low) derived in [1], which give rise to the reentrant be- 
havior, appear as a subset of the more general RG equa- 
tions for XY-systems with additonal symmetry-breaking 
[9] or random fields [10], as well as for solid films with 
quenched random impurities [11], also these systems have 
to be reconsidered, which we will postpone to forthcom- 
ing publications. 
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For the further discussion it is useful to decompose the 
Hamiltonian (1) into a spin- wave part Hsw and a vortex 
part Hv Since the phase transition is governed by Hv, 
wc will omit Tisu; completely. In the continuum descrip- 
tion, the vortex part can be rewritten in the form [1] (for 
simplicity, we set the lattice constant equal to unity) 



/E 
dVQ(r)ln|r-ri| - j^mj. 



(2) 



The integer vortex charges nii satisfy J^i mi = 0. Q(r) 
is a quenched random charge field, which is related to the 
phase shift A(r) by 27r(3(r) = —dxAy + dyA^. Here we 
made the replacement Aij A(r) by going over to the 
continuum description. Since 

[A]d = 0. [A^ir)A(3{r')]d = ad^o (r - r'), (3) 

where [...]d denotes the disorder average, the random 
charges arc anticorrelated. 

The main result of the work of RSN [1] arc the RG-flow 
equations (4a-c) (see also [2], [9], [10]), which describe 
the change of J, a and the vortex number density y al- 
ter eliminating vortex degrees of freedom up to a length 
scale e' 



J2 



dJ 
da 



(4a) 
(4b) 
(4c) 



Here we use the convention, that only the exchange con- 
stant J is renormalizcd and the temperature plays merely 
the role of an unrenormalized parameter. For a = the 
equations (4a), (4b) behave well defined for T — > 0. This 
becomes more clear, if wc rewrite the vortex fugacity as 
y = 6"^"=/-^, where is the (core) free energy of a sin- 
gle vortex on the scale e'. Then (4b) takes the form 
dFJdl = {ttJ ~2T- TT^cr) = 2(r+ - T){T - T_)/T. 

For (7 > 0, the last term on the r.h.s. of (4b) blows up 
at low T, leading to the reentrance transition mentioned 
above. Whereas for high temperatures the 1 /T coefficient 
of the cr term is plausible, since thermal fluctuations wipe 
out the random potential, we do not see a reason that this 
effect could lead to an unlimited growth of the effective 
disorder strength at very low temperatures. Clearly, (4b) 
cannot be valid at zero temperature. Contrary to RSN 
[1], we argue, that the equations (4a), (4b) are valid only 
for sufficiently high temperatures T > T*{a) > T_{a). 

An indication for T* follows from the flow of the vortex 
entropy S'c = -dFjdT, dSJdl = 2-7T^a+Tr{-§^){l- 
2cry). Since dJ/dT < 0, the entropy is reduced for 
T <T* = 2Ja, a < 7r/8, if one goes over to larger length 



scales. This leads finally to a negative entropy, which 
we consider as an artefact of the calculation [1] (see also 
[2], [9], [10], [11]). The vanishing of the entropy in dis- 
orderd systems usually signals a freezing of the system 
by approaching T* from high temperatures [12]. Sim- 
ilarly, the flow of the vortex energy = F^ + TSc, 
dEc/dl = (1 -2(T^)7r( J -T|^) leads for T < T* eventu- 
ally to negative values of the core energy. Inevitably, this 
favours multiple occupancy of vortex positions. However, 
the resulting vortices of higher vorticity jmj > 1 appear 
even in the presence of disorder much less likely than 
those with |m| = 1: since their energy cost scales as 
whereas their energy gain scales only as to. This effective 
repulsion of vortices leads for T < T* to a much smaller 
vortex density than in the RSN-theory [1], which neglects 
completely the interaction between vortex dipoles. 

For T < T* we expect the physics to be different from 
that described by Eqs. (4). Since T*{a) intersects the 
RNS phase boundary at a = tt/S where = = 
j7r/4, the whole (T, (7)-range in which reentrance was 
observed belongs to the freezing region, which has to be 
reconsidered. 

To find the correct behaviour at low temperatures, we 
consider first the system at T = 0. A simple estimate 

shows, that then vortices will not be relevant if the dis- 
order is weak. Indeed, the elastic energy of an isolated 
vortex of charge ±to in a system of radius R is m^wJ In R, 
which has to be compared with the possible energy gain 
Er from the interaction of the vortex with the disorder. 
If we rewrite the second term in (2) as TOjV(rj), we 
find [^^(ri)]^ ~ 2'Kj^a\nR. Hence the typical energy 
gain is — J(TO^27rcrlni?)^/^. 

In order to find the maximal energy gain, we have 
to estimate the number n(R) of vortex positions in 
which the energies VijCi) are essentially uncorrelated. 
Two vortex positions , rj have independent energies if 
[y^(ri)]d >• [V{ri)V{rj)](i , a condition which can be 
rewritten with 

[(y(r,) - V{rj)f]d = 4naJ^ In |r, - r^j = A\r, - r^) 

(5) 

as Injrj - r^j « (1 - e)lni? with e < 1. Thus 
n{R) w i?^^ and the maximal energy gain from ex- 
ploiting the tail of the gaussian distribution for V{ri) is 
Er ~ — 2J(TO^7rea)^/^ Ini?. The total vortex free energy 
at T = is therefore 



« ^7r(m^ - 2(TO2eo-/7r)i/2) Ini? 



(6) 



and hence vortices should be irrelevant for weak disorder 
cr < 1. 

In studying the behavior for T = but larger a we 
have to take into account the screening of the vortex and 
quenched random charges by other vortex pairs. This 
can be done most easily by using the dielectric formal- 
ism. Here we follow the treatment of Halperin [13] who 
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showed, that screening by vortex pairs with separation 
between R and R + dR changes the coupUng constant 
J(i?) (which corresponds to the inverse dielectric con- 
stant) as 



J(i? + dR) = J(i?) - 47r2 ^ a^{R)j'^{R)2TiRp^{R)dR. 

(7) 



m>0 



Here am{R) is the polarizabihty and pm{R) is the proba- 
bihty density of a pair with charge -l-m at ri and charge 
-m at r2, i? = |ri - r2|. 

For the calculation of pm{R) and am{R) we use the 
fact, that the interaction energy between a vortex pair 
and the disorder is gaussian distributed with a width 
A(i?). If the density of pairs is sufficiently small, we 
may neglect the interaction between pairs and write 



-2|m|(Tr,/lnfl+_Ee) 



/2^A(i?) 



27r2m2 Ini? 



R- — 



(8) 



where the r.h.s. of (8) is valid only for ct, Ed J <C ln_R. 
Since p|m|>i(^) ^ Pi{R) =: we will neglect double 

occupancy of vortex positions. 

The polarizabihty a(R) = ai{R) can be calculated in 
a similar way and is found to be a{R) « R? /T* at large 
R. With ^ R'^p{R) and / = In i? we get from (7) and 
(8) for T = 



dj Q J n 

dl a 



dl 
da 

'dl 



= (2 
= 0, 



4ct 



)2/ 



(9a) 
(9b) 
(9c) 



where we again neglected terms of the order a/l. These 
are the flow equations, which replace (4) at zero tempera- 
ture. Within this approximation, the system undergoes a 
phase transition at CTc = 7r/8 from a KT to a disordered 
phase, which is in qualitative agreement with our esti- 
mate (6). At (Tc the exponent rj shows a universal jump 
from 1/16 to zero. For a > ac the y reaches a value of 
order magnitude unity on the scale i? w ^ with 



^ (X e 



l/6(l-7r/8(T) 



(10) 



& is a constant, which depends on the details of the sys- 
tem. For R > ^ our flow equations are no longer valid, 
since y is no longer small. We identify ^ with the corre- 
lation length in the disordered phase. 

We discuss now the properties of the system at low but 
finite temperatures. The T-correction to our free energy 
estimate (6) are of the order — 2rini? (or smaller) and 



hence will not allow a reentrance transition. A more effi- 
cient way for thermal fluctuations to influence the low-T 
behaviour would be the generation of uncorrelated frozen 
charges Q{r). However unlike to random field systems, 
where uncorrelated random fields are indeed generated 
from anticorrelated random fields [14] , which destroy the 
ordered phase in 2 dimension at all non-zero T, we do not 
see such a mechanism here. The main difference consists 
in the existence of a double degenerated ground state in 
the random field system at T = 0. 

The physics at finite temperature can also be captured 
within the dielectric formalism. Neglecting again the 
interaction between vortex pairs at different positions, 
we calculate the normalized probability for a pair with 
charges ±rn as 



PmiR) = P-m{R) 



-E^{R)/T 



(11) 



where E„,{R) = 2m'^{Ec + t:J \nR) + m{V{vi) ~ V{y2)) 
denotes the pair energy. At large R holds p\m\>i{R) ~ Oi 
since the elastic energy cost cx will be compensated 
with decreasing probability by an energy gain oc m due 
to disorder. We therefore drop occupancies \m\ > 1. Fur- 
thermore, for a given configuration of disorder, one of the 
two energies E±i{R) is always so large, that the corre- 
sponding weight factor e^±i(^)/^ can be neglected. After 
this approximation, the probability for a single pair reads 



p{R)=Pi{R) = 



1 



I ^ (.E^{B.}/T 



(12) 



Eq. (|11|) thus effectively reduces to the disorder average 
oi the Fermi distribution function. In other words: vortex 
pairs of vorticity one can be treated as non-interacting 
fermions. In the limit T = 0, where this distribution 
function becomes step-like, Eq. ( [l2| ) immediately reduces 
to the previous expression (8). At finite temperature, the 
disorder average in Eq. (12) is performed by splitting the 
integral over the disorder distribution into two contribu- 
tions corresponding to i?i (i?) ^ 0. To leading order in i?, 
we find p{R) - R-^/'^" for < T < T* = 2J(t, whereas 
p{R) - i?-27r,7/T(i-aj/T) f^j. ^ > T* . Plugging these re- 
sults into the definition of y, we obtain the fiow equation 
(4b) in the whole range T > T*, whereas Eq. (9b) is 
valid in the whole range < T < T* . Both equations 
coincide at the boundary T = T*. 

We add a few remarks: As long as Ei{R) ^ T, the 
Fermi distribution can be replaced by the Boltzmann 
distribution, as is usually done in the treatment of the 
KT transition [4]. The disorder average of the latter 
yields p{R) - ji-^nJ/Tii-aJ/T) Yience Eq. (4b) for 
all temperatures. However, for y < y* tJ^g condition 
Ei{R) 3> T is no longer fulfilled for most of the vor- 
tex positions (see also our remarks below Eqs. (4)) and 
hence this approximation breaks down. Indeed, use of the 
Boltzmann distribution at low temperatures would lead 
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to p{R) ^ 1, and the interaction between vortex pairs 
could no longer be neglected. It is therefore important 
to calculate p{R) from An attempt to improve upon 
the Boltzmann-approximation consists in expanding jl^ ) 
into a power series in e^^^^^-*/^. The n-th order term in 
the expansion yields a contribution /j;-27r,//T(n-n o-j/t) 
to p{R). The series is divergent, i.e. for large R higher 
order terms are more important than lower order terms, 
irrespective of temperature. These higher order terms 
generate contributions to the flow equation dy/dl, which 
tend to blow up y even faster. One might hence expect 
an instability of the ordered phase, similarly to the ob- 
servation of Korshunov [5]. In fact, the above expansion 
and in particular the replacement of the Fermi- by the 
Boltzmann distribution are disqualified a posteriori. 

We conclude, that dy{l)/dl < for aU T < T+. 
The polarizability at finite temperatures is given by 
a = R^/{T + T*) for T < T* and by a = R^/{2T) 
for T >T*. Thus dJ/dl < holds for all T < T+ which 
is sufficient to guarantee the absence of reentrant phase 
topology. In the special case of -Ec ~> oo the phase bound- 
ary is given by T+(o') for T > Jtt/A and a horizontal line 
tJc = 7r/8 for smaller T, as shown bye bold lines in Fig. 
1. This is consistent with the prediction of Ozeki and 
Nishimori [8] about the existence of a horizontal phase 
boundary. We expect the critical behavior at T+{a) as 
discussed in [1] to be unchanged. At finite core energies, 
the actual transition temperature will be renormalized to 
Tc(cr) < r+(cr). Its value Tc{0) is given by the KT flow 
equations without disorder and lies only slightly below 
T+(0) for large Ec- For small a, the critical RG tra- 
jectory flows completely in the domain of equation (4b), 
where weak disorder induces weak additional screening. 
Therefore T'c(o') will smoothly decrease with increasing 
a. We expect this function to end up in Tc{tt/8) = 
monotonously, since flow equations vary monotonously 
in parameter space. 

Our conclusions about the absence of reentrance are 
confirmed also by a discretized Migdal-Kadanoff renor- 
malization group (MKRG) scheme [15] for model (1), 
which we consider in the last part of this paper. Our 
technique has been shown to be similar to that of Jose 
et al. [16]. Their approach is based on studying Migdal- 
Kadanoff recursion relations for the Fourier components 
of the (spatially uniform) potential. 

In the discretized scheme [15] instead of allowing to 
be a continuous variable, we constrain it to take one of 
many discrete values which are uniformly distributed be- 
tween and 27r. Hamiltonian (1) is now defined for values 
of (p restricted to 2TTk/q, where fc = 0, 1, 2, . . . , {q— 1) and 
g is a number of clock states. We define 

Jij{q,k) — J cos{2nk/ q — Aij). (13) 

The recursion relations for Jij{q, k) may be found in [15]. 
For the random 2D system, the numerical procedure is 



based on creating first a pool of Np bonds, each decom- 
posed into q components according to Eq. (p^. One 
then picks Np random batches of 4 such bonds (the cor- 
responding rescaling factor is equal to 2) from the pool 
to generate a new pool of the coupling variables and 
the whole procedure is iterated. We consider typically 
iVp=2000 and g=100. The results depend on these pa- 
rameters rather weakly. 

It should be noted that Gingras and Sorcnsen [16] have 
tried to construct the phase diagram of the 2D random 
Dzyaloshinskii-Moriya model (this model is believed to 
be equivalent to (1) by the same discretized MKRG ap- 
proach. In order to locate the Kosterlitz-Thouless phase, 
they study the scaling behavior of the absolute average 
height of the potential, h, which is defined as follows 

h = (|J,,(q,0)- J,,(g,<?/4)|) . (14) 

Due to errratic behavior of h they could not draw the 
phase diagram. The reason here is that h representing 
only two clock states cannot correctly describe the sys- 
tem with many clock states. 

To obtain the phase diagram one can consider the scal- 
ing of the maximal and minimal couplings for each effec- 
tive bond or the scaling of the average of absolute val- 
ues of all q couplings. The scaling properties of these 
three quantities has been found to be essentially the 
same, so it is sufficient to focus on the maximum cou- 
pling Jmax{q,k). The details of this approach can be 
found in Rcf. [15]. 

It should be noted that the discretized MKRG ap- 
proach cannot rigorously reproduce the quasi-long-range 
XY order in 2D [15]. The scale invariance of Jmax{q,k) 
in the KT-phase is merely approximate in this approach. 
In practice, the scale invariance of Jmaxil, k) persists for 
about 20 iterations. Further iterations lead to an even- 
tual decrease of JmaxiQ^k) at any nonzero T. Having 
this caveat in mind, we can locate the boundary between 
the paramagnetic and KT- phase (see Fig. 2). Thus the 
MKRG gives us additional evidence that the reentrance 
is absent in model (1). 

To conclude, in the present paper we have shown by 
a combination of simple analytical arguments, a renor- 
malization group calculation and a Migdal-Kadanoff RG 
scheme at finite T, that the 2-dimensional XY-model 
with random phase shifts does not exhibit a reentrant 
transition. 
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FIG. 1. (ct,T) phase diagram of the model (1). T±(cr) 
are the upper bounds for the transition temperatures Tc(a) 
and Tre(o-) between the disordered and the KT phase in the 
RSN-theory [1]. Note, that r_-line lies completely in the 
freezing region (hatched area). The true phase transition line 
Tc(cr) is denoted by the dashed line which is bounded by T+{a) 
and a = -k/8. The line Tre is not shown here. 

FIG. 2. (cr, T) phase diagram obtained by the discretized 

Migdal-Kadanoff RG scheme. PM and KT denote the para- 
magnetic and Kosterlitz- Thouless phase respectively. In the 
PM phase Jmax (?, k) scales down monotonously whereas in 
the KT region it reaches a fixed value at large scales. The 
critical values ^^^(T = 0) « ksTda = 0)/J « 0.44. One 
can also demonstrate that the phase diagram of the random 
2D Dzyaloshinskii-Moriya model has the same topology. 
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